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(N ■ Abstract 

^t- : 

This paper shows that a hyperbolic equation for heat conduction can be obtained directly 
using the tenets of linear irreversible thermodynamics in the context of the five dimensional 
space-time metric originally proposed by T. Kaluza back in 1922. The associated speed 

O 

of propagation is slightly lower than the speed of light by a factor inversely proportional 

U '. 

to the specific charge of the fluid element. Moreover, consistency with the second law of 

£> ■ 

• — h | thermodynamics is achieved. Possible implications in the context of physics of clusters of 

' 

$— i ' galaxies of this result are briefly discussed. 
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I. INTRODUCTION 



It is today a well established fact that knowledge of magnetohydrodynamics 
(MHD) is essential in our understanding of a vast number situations which occur 
in astrophysics l|. Less known is the relationship which exists between the basic 
equations of MHD with non-equilibrium thermodynamics, specially in the case in 
which one wishes to include dissipative effects into the scheme. Four years ago, 
two of us wrote a paper showing clearly that when this is the case, the classical 
non-relativistic structure of linear irreversible thermodynamics (LIT) restricts the 
constitutive equations relating the electrical current to an electric field to the well- 
known Ohm's equation when the second law is satisfied. Non-ohmic effects can thus 
be included by attempting an extension of LIT to include such cases, but a theory 
with these characteristics which is also at grips with the second law of thermody- 
namics does not yet exists ^|. The other alternative is to approach the problem 
by incorporating the effects of the electromagnetic field using a relativistic multi- 
dimensional theory, a step that also brings the theory to the framework of general 
relativity. Many multidimensional theories have been formulated since Einstein's 
idea of unifying nature's fundamental forces in a single theory was put forward. One 
of the pioneers in this topic was the mathematician T. Kaluza Although he was 
not the first one to propose an extended metric, his idea of adding extra elements to 
the 4x4 general relativity space-time can be considered as the origin of a variety 
of multidimensional theories that has been proposed in the last five decades. In this 
work we formulate the thermo-hydrodynamics for an ionized fluid in the context of 
Kaluza's original five-dimensional theory. This framework, which does not account 
for weak and strong interactions, has sufficient information to include gravitation 
and electromagnetism. 

In Re, Q it was shown how the ~e a risi n g ,o m the extra dimension a. 
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counts for the electromagnetic effects in the hydrodynamic equations. This was 
achieved by using Meixner's theory of non-equilibrium thermodynamics j^, Q] and 
the 5x5 metric proposed by Kaluza. In that work the entropy production was 
calculated and it was shown how, using this formalism, the constitutive equation for 
the electric current can be other than Ohm's law, which as we know is not always 
valid. However, the causality issue in the heat equation was not addressed. 

In this work we complete the calculation in Ref. P| and add a new fundamental 
element in Kaluza's magnetohydrodynamics. A parameter £ is introduced in the 
metric which ultimately adds non-negligible corrections in the transport equations. 
The causality issue in the transport equations is also addressed by calculating the 
heat equation, which is shown to be of the hyperbolic type, and thus does not allow 
heat waves propagating with a speed larger than the speed of light. 

Section Ull is a brief summary of Kaluza's formalism. In Sect. lIIII we establish the 
entropy production and the heat equation showing how causality is not violated if the 
five dimensional formalism is considered. This is accomplished by using the standard 
projector h£ || which eliminates the somewhat delicate concept of dissipation in 
time A summary and a brief discussion of the results here obtained is included 
in Sect. lIVl In Appendix A some relevant details of Kaluza's theory are shown 
including the relevance of the parameter £. 



II. KALUZA'S FORMALISM 

The theory formulated by T. Kaluza 0|, unifies general relativity and electromag- 
netism by extending the 4-dimensional manifold into a 5-dimensional space-time. 
In this metric, the fifth column contains the components of the electromagnetic 
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potentials: 
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H 

r. ^> 



(1) 



where are the components of the vector potential and 4> is the electrostatic scalar 
potential. The constant £ is a parameter whose value is set equal to 



16irGe 



(2) 



by matching Einstein's field equation to a Poisson equation for 0, in the non- 
relativistic limit, as is shown in Appendix A. For simplicity we will neglect the 
magnetic field and assume that the 4-manifold is given by a quasi-Minkowski metric. 
That is, if the gravitational potential is ip, the metric tensor can be written as 



and 
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neglecting terms of order 1/c 3 and lower. In this framework, the Christoffel symbols 
are proportional to the components of the forces and 

dx 5 q 

dt ~ mC [ ) 

so that the equation of motion for a charged particle in a gravitational and electric 
field corresponds to a geodesic in this space-time as shown in Appendix A. It is 
also required that quantities have no variation in the fifth dimension by setting 
d/dx 5 = 0. This is called the "cylindrical condition" Q and follows the idea that 
the extra dimension is closed on itself as is now assumed in some multi-dimensional 
theories. 



III. HEAT CONDUCTION IN KALUZA'S MAGNETOHYDRODYNAMICS 

As shown in Ref. P|, the MHD equations can be obtained within the five- 
dimensional theory. The effects of the electrical force arise in a natural way from 
the curvature of space-time. No external forces are introduced since the components 
of the electromagnetic field are included in the equations through the Christoffel 
symbols in the covariant derivatives (see Appendix A). 

The conservation of particles is unaffected by the fifth dimension. That is, if n is 
the local particle density 

n + n0 = 0, (6) 

where 9 = is the divergence of the velocity for which the last term, w 5 5 , vanishes 
due to the cylindrical condition. For simplicity, we will assume that the calculations 
are performed in the comoving frame. Thus, 6 = and the particle conservation law 
is given by h = 0. 
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The stress tensor within this framework is given by P| 

T^pu^+pK + E*, (7) 

where p is the mass-energy density, p the hydrostatic pressure and H£ is the viscous 
stress tensor. The spatial projector h$ is defined in the standard way so that 

u v K = 0, (8) 



and, as will be seen, 



This is achieved by having 



u^ v = 0. (9) 



K = 6j£+ — 2 u v u\ (10) 
or 

where a 2 = c 2 —ugii . Some authors include a heat term in the definition of the stress 
tensor following Eckart's formalism In Meixner's theory, heat is viewed as a non- 
mechanical form of energy and is included in a separate equation. A discussion on 
the difference between both formalisms and a possible way of discerning the correct 
way to treat heat in general relativity via a light scattering experiment can be found 
in Refs. 0, In Meixner's thermodynamics, one considers the total energy flux: 

Jfa = u v T* + n e int u» + Jf Q] , (11) 

where u v T^ and ei nt u^ are the mechanical and internal energy fluxes respectively 
and Jjgj is the heat flux. Total energy conservation requires Jm.^ = which, using 
the form of the stress tensor given by Eq. 0, yields 

= «*ss. (12) 

To obtain Eq. (|T2l we used that, according to Einstein's field equations, = and 
that u^u^ = a 2 is a constant. Thus, in the comoving frame (m^)^ = 0. Hence, the 



balance equation for the internal energy is 

ne int ^-u^-J^. (13) 

The electric contribution to the internal energy is, as mentioned above, contained in 
the covariant derivatives of Eq. (fl~3|) . It is pertinent to point out that Eq. (fl~3|) would 
differ from its Eckart 's type counterpart by a term proportional to the hydrodynam- 
ical acceleration. In both cases, in the comoving frame, both formalisms would lead 
to Eq. (EH). 



A. Entropy production 



In Meixner's theorvthe entropy production is obtained by using the local equi- 
librium assumption |5j, |6j namely 



s = s{n, e int ) , (14) 
for which, since h — 0, one can write 



n 

ns = —e int , (15) 



or, introducing Eq. {T 



ni+ (¥) =4«r,s;-|lr,„. (16) 

One can identify the right side of Eq. (jTfil) with Clausius' uncompensated heat or 
the entropy production which, according to the second law of thermodynamics is 
required to be non-negative: 

J[Q] -T, - u»*>t > 0. (17) 



J^2 
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This requirement suggests constitutive equations as expressions of each source term in 
Eq. (fT6|) as a product of a thermodynamic force and its corresponding flux. Following 
this reasoning we write a Fourier law for the heat flux: 



Jf Q] = -Khp*. (18) 

In Eckart's theory a second term, proportional to the acceleration, is included in 
Eq. (JT8l) which is not present here. This is not only because u % = in the comoving 
frame but primarily since, as mentioned above, the heat flux is treated separately in 
Meixner's formalism [ijj]. Anyway, both formalisms yield the same expressions for 
the entropy production, heat flux and transport equations in the comoving frame. 
In order to guarantee that the second term in Eq. 1)17)) is positive we write 

^ = -fltf, (19) 

where r\ is a "viscosity" coefficient. Equation ()T9|) relates the viscous stress tensor with 
the symmetrical part of the (traceless) velocity gradient a£. We want to remark that 
Eq. (P2J) in the comoving frame is equivalent to the standard constitutive equation 
E£ = —rjh^h^a^. Bulk viscosity has been neglected. A calculation for in the 
comoving frame leads to (see Appendix A): 



for i — 1, 2, 3 and the rest of the coefficients vanish. Since in Einstein's equation, 
G^is symmetrical, the viscous stress tensor does not have antisymmetrical compo- 
nents and thus we can assume that the coupling constant which would correspond 
to a rotational viscosity vanishes. Then, 

E) = for i, j = 1, 2, 3, 5, (21) 

and 

^ = + (22) 



When these constitutive equations are substituted in Eq. (fT6|) . the entropy balance 
equation reads 

» + (f) ^ = + ,± (*, + (r + ■ (23) 

The second term in the right side of Eq. is due to the fields present in the problem 
and is obtained here as a consequence of the curvature of space-time. 



B. Heat equation 

The heat equation is obtained by using the local equilibrium assumption for the 
internal energy 

dT 

ne int = nC v —, (24) 

together with Eq. (fT3|) and the corresponding constitutive equations given by 
Eqs. (fT9|) and p2|) . The first three terms in the divergence of the heat flux give 
the standard result 

(jy ;4 = - K V 2 T fori = 1,2, 3. (25) 

Note that, in the standard theory, the heat flux only contains the three terms in 
Eq. (|2*51) which, when introduced in Eq. (1231) . yields a parabolic heat equation. How- 
ever, in the five-dimensional theory, the fourth component of the heat flux is nonzero 
retaining its consistency with the orthogonality condition in Eq. (jHJ. That is, since 

hf = for i = 1, 2, 3 (26) 



hi = 1 - (27) 

and by means of the cylindrical condition T' 5 = 0, the fourth component of the heat 
flux is 

4n = -« (i - t4^) t ' 4 > (28) 



where 5 2 = 16 ^ e Taking the covariant derivative in Eq. (f28|) leads to 

The nonlinear term • VT couples the gravitational field to the temperature gra- 
dient. It is worth mentioning that there is also a coupling of order 1/c 4 between 
the electrostatic field and the temperature gradient which consistently has been ne- 
glected. An analysis of the consequences of these couplings is an interesting topic 
which will be addressed in the future. With these extra terms, Eq. can be written 
as 

*§ - ^4 -rh)(w- *rf (* + **) + ^) 

(30) 

where D T = nC v /n is the thermal diffusivity. Equation (|30|) is a hyperbolic, and 
thus causal, heat equation. Since the coefficient rj plays the role of a viscosity, the 
factor that multiplies the dissipative source on the third term on the right side can 
be interpreted as a Prandtl number that quantifies the ratio between viscous and 
thermal dissipations. Diffusion dominates for small values of this parameter and, 
assuming the fluid is neutral, the standard heat equation is recovered. Otherwise, 
since e$G is a very small quantity, even a very small net specific charge will yield 
5 2 > 1. That is, because of the large value of l/£ the corrections due to the electric 
charge of matter, and thus to the fifth dimension, are always significant unless q/m 
is identically zero. 

The speed of propagation for the wave-like solution of Eq. (JH5j) is lower than the 
speed of light. If c is the modified speed and since 5 2 > 1 we obtain 




As an example, consider an element of fluid purely constituted by electrons. In that 
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case, the specific charge is e/m e and thus 

c 2 = c 2 (1 - lO" 42 ) . (32) 

Having 5 2 ^ is the key element for having causal heat conduction, although the 
correction imposed on the speed of propagation is small. Therefore, this will not 
make significant changes in the measurements where one would obtain c ~ c. 



IV. SUMMARY AND DISCUSSION 



The main objective of this paper is to show how using the unifying ideas set forth 
by Kaluza over eighty years ago and the tenets of LIT one can derive transport equa- 
tions in MHD which incorporate not only electrical dissipation but other dissipative 
effects. This is clearly exhibited in Eq. pUj). the main result of this paper in which 
the heat equation is modified in two ways. Firstly, it is a hyperbolic equation which 
therefore does not violate causality and is at grips with the second law of thermody- 
namics. Secondly, it includes coupling effects between the external and temperature 
fields together with viscosity. Thus, a relativistic framework allows for the inclusion 
of dissipation consistently with a non negative entropy production. Clearly, in the 
limit when c — > oo one recovers the standard equation for heat conduction. 

Finally, it is worth mentioning that the usefulness of these results is not merely 
conceptual. There is an increasing interest in dealing with the inclusion of dissipative 
effects in the physics of intracluster media [ll| and/or clusters of galaxies jl^] when 
temperatures are close to 10 8 K. The existence of significant magnetic fields in the 
intracluster medium is a fact that encourages further applications of this formalism 
of relativistic MHD. 

This work has been supported by CONACyT project 41081-F and FICSAC, Mex- 
ico (PFSA 2005-2006). 

11 



Appendix A: KALUZA'S THEORY AND THE PARAMETER f 



The Christoffel symbols are defined as 



r$ y = 9 — ^ + ^-^\. (Ai) 



Pl 2 \dxn dx? dxv 

Using the metric tensor defined in Eqs. ((H) and (jHJ), with Ai = one can easily obtain 
that 



ri 5 = ±<f>,i n 4 = (A2) 



r* - 1 
2c^ i44 " c 2 

The "shortest path" in this five dimensional space corresponds to a geodesic which 
is given by 

g + c 2 r>V = 0. (A3) 

It is easy to show that using the elements in Eq. 1(A2J) and u 5 = q/m£, Eq. (|A3jl 
corresponds precisely to the equation of motion for a charged particle that moves in 
a region where the gravitational and electrostatic potentials are ip and <fi respectively. 
Also notice that, since the covariant derivative is defined by 



dV h 

;M dx^ w y ' 



for any tensor V v , the components of the gradient of the potentials will appear in a 
natural way in the equations. For example, when computing a given component of 
the velocity gradient in the comoving frame 



u. 



r>" = r> 4 + r\ 5 u 5 , (A5) 



which, when the relations given in Eq. (|A2jl are introduced, yields Eq. (|20|) . 

The electromagnetism equations in this framework can be obtained from the field 
equations. That is 

-R44 - -^Rgu = kT 44 , (A6) 
12 



yields the Poisson equation for the gravitational potential with k = while 

i?45 - 7^45 = KT 45 , (A7) 

corresponds a Poisson equation for the electrostatic potential. As usual, R^ u repre- 
sents the Ricci tensor. The coefficient £ is fixed by comparing this equation with 

V 2 0=--- (A8) 

This part of the calculation is subtle and deserves a closer look. In Kaluza's original 
work, the curvature scalar R in Eq. (|A7j) is neglected. However, calculating i? 45 and 
using the relation 

r = -«r;, (A9) 

one obtains, for the field equation, 

|r* + H* -;£#=-?*» < Aio > 

where p e = pq/m is the electric charge density. In this Helmholtz-type equation, the 
second term corresponds to the scalar curvature term in Einstein's field equation. 
Rearranging terms, one can compare this linear term with the source on the right 
side as follows. By assuming that the second term can be neglected one can obtain 

/I6p^ 6 x 10 -i9 % .A. (A11) 
V cr cm 

Introducing this estimate for £ in Eq. IjAlOjl yields 

W ~ U± - 2c) f p., (A12) 

where we neglected the first term in the parenthesis in Eq. (|A10|) since c 2 <C ~ 
^2 x 3 x 10 36 (-^ • ~) . Under the assumption that <f>/c is small, one can also neglect 
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the linear term in and the standard Poisson equation for the electrostatic potential 
is obtained. 
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